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Abstract 

A new exponent characterizing the rounding of crystal facets is found by mapping a crystal surface onto 
the asymmetric six-vertex model (i.e. with external fields h and v) and using the Bethe Ansatz to obtain 
appropriate expansions of the free energy close to criticality. Leading order exponents in 5h, 8v are determined 
along the whole phase boundary and in an arbitrary direction. A possible experimental verification of this 
result is discussed. 
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Along a first-order coexistence line the shape of a solid inclusion in equilibrium with its fluid is determined 
by the Gibbs-Curie constrained minimization condition: for a fixed volume V of the inclusion, the equilibrium 
surface Q = z(x,y) is such as to minimize the surface free energy 

F(T) — / r(n, T) dfi = minimum (1) 

where the surface tension r depends on the orientation oiQ.fi — V z(x, y) and temperature T. To solve the 
constrained variational problem to determine z(x, y) one introduces a Lagrange multiplier 2v 

4 [[f(p,q)-2vz(x,y)]dxdy\=0 



P = §x> Q = TFy are ^ ne slopes of the surface and f(p, q, T) — r(p, q, T)y/l + p 2 + q 2 is the specific surface free 
energy projected onto the (x, y)-plane. The solution reads fiL || 

vz = f(-vx,-vy) (2) 

where the Legendre transformed potential / is defined through (|] 

7{Vx,Vv) = f(p, l) ~ PVx - qv v (3) 

with r) x — yt, T) y — -7jt the surface tilting fields conjugate to the slope variables p and q. Thus the determi- 
nation of / yields the equilibrium surface z(x, y) directly, up to some overall irrelevant constant. The main 
question however is how much of the physics can be captured by this description. Crystals are believed to 
be entirely faceted at T = 0. As the temperature is raised, thermal fluctuations bring about a rounding 
of the edges where planes with different indices (hkl) meet. At the roughening transition temperature T R 
the remaining planar portion of (hkl) disappears entirely and the face becomes rough. Different facets may 
have different T R 's. It is clear that questions concerning the exact value of T R or exponents characterizing the 
macroscopic smoothing out of the edges cannot be addressed in a satisfactory way unless one has a microscopic 
model from which r can be derived from first principles. This problem has been tackled with success by several 
authors mainly by using a connection between the theory of equilibrium shapes and two-dimensional lattice 
models, with the asymmetric six-vertex model playing a major role S ||, |sl q|. This system of interacting 
dipoles on a lattice is the natural extension of the symmetric problem through the inclusion of external fields 
h and u 0. It is a classical paradigm of models soluble through the Bethe Ansatz and an exact expression 
for its free energy has been known for many years (qj. Originally introduced in the context of equilibrium 
crystal shape problems (ECS) to describe the (001) meet of a bec crystal 0, it was later extended to cover 
other facets on different lattices as well. In this way very comprehensive results have been obtained regarding 
the form of the equilibrium surface as well as its thermal evolution |5| . In particular there is a region in the 
(h, v) plane bounded by a curve T and containing, the origin (h = 0,v = 0), which defines the loci of points 
where the free-energy retains its zero-field value [|S| . This implies that z is a constant over an entire region of 
the parameter space. Beyond V the free energy changes continuously as a function of the fields. Appropriate 
expansions of the surface free energy around T show a leading order exponent for the scaling of z(x, y) as one 
approaches the boundary of the (110) fee facet along the y = direction |0 



z~ {x-x c )* +0{x-x c y (4) 

This result has been confirmed experimentally M, [Io|]. However, even though a great deal has been learned 
from exact analytical methods, their handling is no trivial matter and previous results were somewhat limited 
in their scope by rather involved technical difficulties. In spite of the fact that it had already been conjectured 
that (ji|) should hold on the entire facet boundary a proof of this result still eluded us ||. In this letter 
we give such a proof based on the expansion of the free energy of the asymmetric six-vertex model via the 
Bethe Ansatz, introduced in 11, [l2| as a generalization of previous methods §. The main results can be 



summarized as follows: the exponent | dominates the rounding in all directions but the tangential one, 
which case a new exponent 3 dominates. These results hold for the whole boundary. 

For the sake of clarity, we consider one particular geometry, namely that of the (110) facet of an fee 
lattice |lj| . One may associate an energy to the surface arising from the cleavage of an infinite ideal crystal 
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through the plane (hkl) based solely on the number and type of bonds broken per primitive cell. This is the 
simplest model possible and yet one which allows an exact handling of the equations that follow. Consider 
an fee lattice (fig. 1) with energies —J\ and — Ji assigned to nearest- and next-to-nearest neighbor bonds 
respectively (J\ > J 2 > 0). The T = surface tension reads Q 

r(n hkl ) = *™ 1 • ( Ji[210] + MU1}) (5) 
V p \nhki\ \ J 

where rihki = hx + jy + kz is the vector normal to the plane [hkl) and V pc is the volume of the primitive fee 
cell. Voids (excitations in the bulk) as well as overhangs (height differences between neighboring columns of 
atoms bigger than some unit length) are strictly forbidden (solid-on-solid condition) . The connection between 
the (110) facet and the vertex model is shown in figure 1. From (||) one may associate the following energies 
to the vertices || 

ei = J2 + - =J 2 -H-V 

Vx 

e 2 = J 2 = -h + H + V 

Vx 

e 3 = Ji —riy = Ji - H + V 

v2 

e 4 = J 1 H — -j=i] y = Ji + H — V 
v2 

e-b = e 6 = (6) 
Here A = ^ is the area of a vertex (for a lattice constant a) and 
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Note that this parametrization is slightly different than the one used in [|| but is consistent with figure 1. 
The free energy (|J) follows straightforwardly Q 

J(ff, T) = /(HO) + f 6 - vert ex(Ji, J 2 ,H, V, T) (8) 

A configuration of alternating vertices 5 and 6 corresponds to the macroscopically smooth (110) plane. The 
twofold degeneracy (exchanging vertices 5 <-> 6) represents the invariance of the crystal surface by the removal 
of the top layer of atoms. In terms of the original dipoles, this configuration of vertices has a zero net dipolc 
moment (which trivially translates to the facet (110) having zero slope, being the reference plane). Regions 
of net dipole moment different from zero (nonzero slopes p, q) correspond to tilts away from the reference 
plane. To find / (or equivalently /) of this system of interacting dipoles attached to the links of a square grid, 
we consider a lattice of N x M sites and impose periodic boundary conditions on both directions. Define the 
row-to-row transfer matrix 

N 

^) (2 }m = EI1C; i m ( 9 ) 

{£} fe=i 

where {a} = a%, • • • , an are the arrows on a row of TV vertical links (the a n s are one row above) and [3k 

is the arrow to the left (right) of the k-th vertex. R^ k J^ +1 is the Boltzmann weight associated to a given 

configuration {ak, a' k , (3k, Pk+i} of the k-th vertex (see fig. 2). It is more conveniently described in terms of 
the so-called spectral parameter u and anisotropy 7, whose explicit relation to the original parameters reads 



sinhit = e~ Jl/fcBT \/cosh 2 7- 1 (10) 
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cosh 7 



cosh 



_i_ (Jl _ j2) + ^».+.« 



(ii) 



and k B is the Boltzmann constant. We also incorporate the Boltzmann factors in new definitions for H and 
V, namely h = f3H and v = (3V (see fig. 2). 

The transfer matrix can be diagonalized exactly with the Bethe Ansatz. Its eigenvalues read M 



A(u,7,h,«) 



e v(N-2n) e hN 



sinh(7 - 



+e v(N-2n) e -hN 



sinh7 

sinhit 
sinh 7 



JV n 



■Q sinh (i 



J= i sinh( 2 

N n _•„!./— 37 



n 



sinh(- 



=1 sinh(^ 



(12) 



where the are the roots of the Bethe Ansatz equations 



sinh(^ 



2 > 



sinha _ *a) 



N 



= (-1) 



-1 2hN 



n 



sinh(7 + ±(a k - afi) 
sinh(7 - Ua k - at)) 



k = 1, 2, • • • , n 



(13) 



n represents the number of arrows in a row of N vertical links which are reversed with respect to the reference 
ferroelectric state |f"f • • • f). 

The specific free energy is given by / = — lim^v.M— ><x> ■jarf l n Z where Z = Tr(T M ) is the partition function. 
In the thermodynamic limit this expression is dominated by the largest eigenvalue Ao of T and and the free 
energy reduces to 

j(u,7, = - hm fc R T- 



JV- 



A 



At this point it is interesting to discuss some already known facts about the physics of the asymmetric six- 
vertex model before solving the Bethe Ansatz Equations explicitly. The phase diagram and the nature of the 
phase transitions are well understood when h = v = (symmetric six- vertex), or when h — and u ^ ||. 

The T = phase diagram for (h, v) ^ is trivial In this case the ground state wanders, as a function 
of the external fields, through each one of the different regions composed exclusively of vertices of type 1, 
2, • • •, 5 + 6. The boundaries between each region are sharply defined. In figure 3 a plot of the free energy 
both in the (h, n)-plane as well as in the (t] x , r) y )-pl&ne are presented. The crystal facets to which each region 
corresponds are accordingly indicated || . 

As the temperature is raised, thermal excitations bring about a rounding of edges ||. When T ^ 
(7 finite) the free energy f(h, v) remains constant as function of the field ('flat phase') in a region of the 
(h, v) plane containing h = v = and bounded by a curve T. Beyond this, the field is sufficiently strong 
to destroy the antiferroelectric order of the system, but not strong enough to impose ferroelectric order, and 
an incommensurate phase appears where the polarization (or surface slope), which is zero in the flat phase, 
changes continuously with the field. The spectrum of the transfer matrix is gapless with finite size corrections 
typical of the gaussian model (l5). The curve V has been investigated only at a few points Q| and the 
nature of the phase transition was found to be of a Pokrovskii-Talapov type |p.6[ . To study this curve in 
more detail, one has to solve the Bethe Ansatz Equations in their limiting fornifor arbitrary values of h, v 
and determine the free energy. This is a non trivial problem. A detailed exposition of the methods used to 
tackle |[|) can be found in || O, ||. Here we point out only the main results which are of relevance to the 
questions addressed. Since the Boltzmann factor is inessential to the forecoming discussion, we drop it out of 
our analysis. The vertical polarization y is defined through 



v = lim = lim ( 1 ) 

AT^oo N N^oo\ N J 



(14) 



Inside the boundary curve T in the (h, u)-plane the polarization is constant at its y = value; the free energy 
does not depend on the field and is given by S 



f(u, % h,v) = -2j2 



,-2-yn 



' n cosh 771 



sinh(nii) sinh 71(7 — u) 



(15) 
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The parametric equation (h(b), v(b)) of T is given by |], [TTJ 



M6) (16) 

v ' 2 ^ y ' ncoshn 7 v ' 

n—1 ' 

and 

v( b ) = 7 _ 17 - 2« - 6| g (-T sinh K 7 - | 7 - 2u - ftj)] 
w 2 2 ^ n coshn 7 ' ~ ~ ' v 7 

n— 1 

These equations reproduce only half of the curve T. The other half can be recovered from the symmetry 
f(—h,—v) = f(h,v). Figure 4 depicts T in the (h, t>)-plane along with the same plot in the (j] x , ?]j,)-plane 
(the latter being the boundaries of the (110) facet) obtained numerically for cosh 7 = 21, u = 0.5 (non zero, 
finite temperature). Two points, (h c ,v c ) and (— h c , — v c ) are singled out: the reason for this is explained in 
what follows. 

It is believed ||, |5|, [| that the free energy singularity should be governed by a universal exponent |. A 
calculation performed for i] y — gives || 

/ ~ Ji + constant x (r, x - rf^f' 2 + 0(S V 2 X ) 

From this result, it follows straightforwardly that 

z = zo + constant x (x — xo) 3 ^ 2 

An earlier result also support this observation, the direction of approach being however tilted relative to the 
normal of the surface |Q . It remains however to be proven that such results apply irrespective of the direction 
of approach and for the whole facet boundary. In order to prove this, one has to perform expansions on the 
free energy by means of an extension of previous techniques so as to allow an arbitrary angle of approach 
[p"H |T^| . Consider first the points (h c , v c ). Here the equations simplify significantly and the increments Sh, Sv 
are perpendicular and tangential to the boundary, respectively. One obtains, when approaching T from the 
incommensurate phase 

/ 2 x 3/2 

f(u,j,h c + 8h : v c ) = f(u,"/,h c ,v c )-2c 2 {—8h 

f(u,1,h c ,v c + 5v) = f(u,f,h c ,v c ) - 3( |^ ) |<5u| 3 

where c\ and ci are functions of 7 and u only. These results are not related to the unique geometrical character 
of the points. They hold over the whole boundary curve and can be expressed using the parametric equation 
of r, (h(b),v(b)) given by eqs. (16) and (|l7|). We introduce the following notation fll^ ] 

Mb) - j h v{b) h t {b) = j b h{b) 

vi(b) = j h v t {b) h 1 {b) = j b h t {b) (18) 
A = v t hi - hfVi 

and approach the curve T from the outside along straight lines h — h{b) + Sh, v — v(b) + Sv where Sv = kSh 
and k fixes a slope not tangential to T. The free energy singularity reads in this case |12| 

Sf(Sh) = (^ 2 ^[2(k-^)Sh]V 2 (19) 

These results clearly show that the exponent 3/2 is obtained for all directions except the tangential one, that 
is when Sv = j^Sh. Here the singularity is 



c 
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where c+ 7^ c_ are constants depending on the parameter b of the curve. 

The result of our long calculation is that one has two critical exponents 61 = 3/2 corresponding to the 
rounding off of a crystal facet when the facet boundary is approached in all other directions except the 
tangential one and 62 = 3 if the facet boundary is approached tangentially. Experiments on small Pb-crystals 
[p| which looked for 9\ have found Q\ — 1.60 ±0.15 whereas experiments on large He 4 -crystals pi have found 
ffi = 1.55 ± 0.06 confirming the theoretical predictions. According to A. Babkin |Hj and V. Henrich jL8| 
measurements to put in evidence the second critical exponent O2 are feasible. An open problem common to 
the measurements in both exponents is the effect of gravity. These effects have been discussed in |l9| . 

The authors would like to thank V. Rittenberg for his continuous encouragement and A. Babkin and V. 
Henrich for sharing their expertise on crystal surface experiments. We would also like to thank R. Shankar 
and D. Kim for the fruitful discussions. SRD would like to thank the hospitality of the Physics Department 
at Yale and CAPES (Brazil). 
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Figure Captions 



Figure 1. The fee cell and vertices. On the top left the (110) plane is depicted (full circles). Atoms in the 
plane immediately above are represented by empty circles. Two such planes combined define the 
macroscopically smooth (110) crystal surface. Vertices (full lines, top right) are defined for each group of 4 
atoms (two filled and two empty circles). In the middle, a particular configuration of 4 vertices is depicted, a 
— sign indicating an atomic vacancy. The exact correspondence between vertices and configurations of the 
crystal plane is shown in the bottom figure. Vertices 1 through 4 represent tilts away from the < 110 > 
direction, and correspond locally to planes (111), (HI), (100) and (010) respectively. An alternating 
configuration of vertices 5 and 6 correspond to the (110) surface. 

Figure 2. The vertices and their corresponding Boltzmann weights. (3 is the usual Boltzmann factor {k^T)^ 1 . 

Figure 3. T — view of the free energy in the (h,v)- and (rj x , 77 y )-planc. The values of the parameters are 
A = B = ^±i, C = Ji, D = J 2 and E = \/2{Ji - J2); A is the area of a vertex. 

Figure 4. Parametric plot of T for cosh 7 = 21 and spectral parameter u = ^ in the (h,v)- and 

(Vx, %)-plane respectively. The latter corresponds to the boundaries of the (110) fee facet. The points 

(h c ,v c ) (lower half-plane) and (— h c , — v c ) (upper half-plane) are singled out. 
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